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INTRODUCTION 


In  recent  years,  research  has  been  conducted  in  the  determination 
of  precipitation  Intensity  from  geosynchronous  satellite  Information. 
Studies  by  Sikdar1,  Stout,  Martin  and  Sikdar2,  Schofield  and  Oliver3, 
and  Adler  and  Fenn4  dealt  specifically  with  determining  convective 
rainfall  rates  in  tropical  areas  using  satellite  infrared  (IR) 
brightness  data.  The  majority  of  this  work  has  been  related  to 
inference  of  tropical  convective  precipitation  with  much  less  work 
being  done  to  develop  similar  techniques  for  continuous  precipitation 
areas  in  middle  latitudes. 

Convective  bands  or  clusters  imbedded  within  areas  of  less 
vertically  developed  stratiform  cloudiness  are  commonly  observed  in  the 
cold  sector  of  rapidly  deepening  cyclones  as  well  as  mature  occluded 
storms.  These  imbedded  convective  areas  are  nearly  always  associated 
with  increased  precipitation  intensity  relative  to  areas  of  surrounding 
stratiform  cloudiness. 

In  an  attempt  to  aid  development  of  simple  techniques  applicable 
in  middle  latitudes,  this  study  addresses  the  inference  of 


Ti  Sikdar,  D.  N.,  1971:  ATS-3  observed  cloud  brigtness  field  related 
to  a  meso-to-subsynoptic  scale  rainfall  pattern.  Tellus,  24, 
400-412. 

2.  Stout,  J.,  D.  Martin  and  D.  N.  Sikdar,  1979:  Estimating  rain  with 

geosynchronous  satellite  images.  Mon.  Wea.  Rev.,  107,  585-598. 

3.  Schofield,  R.  A.,  and  V.  J.  Oliver,  1977:  A  scheme  for  estimating 

convective  rainfall  from  satellite  imagery.  NESS  Tech.  Memo.  86, 
NOAA,  47  pp. 

4.  Adler,  R.  F.,  and  0,  0.  Fenn,  1978:  Thunderstorm  intensity  as 

determined  from  satellite  data.  0.  Appl.  Meteor.,  18,  502-517. 
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precipitation  Intensity  and  Its  temporal  change  associated  with 

Imbedded  convective  activity  at  a  location  In  areas  of  continuous 

stratiform  precipitation  from  Geosynchronous  Operational  Environmental 

Satellite  (GOES)  IR  digital  data. 

Imbedded  convective  activity  is  easily  identified  in  high 

resolution  visible  GOES  Imagery  if  high  cloud  layers  are  not  present  to 

obscure  the  convection.  At  low  sun  angles  with  the  shadowing  effect  of 

the  convective  towers,  details  of  these  convective  areas  are  especially 

enhanced.  However,  with  the  need  for  these  methods  to  be  applicable 

during  both  daytime  and  nighttime  hours,  only  the  IR  11  micron  channel 

data  were  used  in  this  study  which  was  based  on  several  synoptic  cases. 

In  each  of  the  cases  examined,  visible  and  IR  imagery  and  loops  were 

used  to  help  pinpoint  optimum  areas  for  study. ^ 

Arrays  of  digital  IR  brightness  values  were  obtained  through  the 

Man-computer  interactive  Data  Access  System  (McIDAS)  of  the  Space 

Science  and  Engineering  Center  at  the  University  of  Wisconsin-Madison.5 

These  arrays,  approximately  36  nautical  miles  (nm)  on  a  side  with  324 

data  points  each,  were  centered  on  surface  observing  stations  which 

experienced  highly  variable  precipitation  intensity  for  a  period  of 

several  hours  during  each  case.  Five  numerical  parameters  were  then 

computed  from  each  array  to  test  their  utility  for  inference  of 

precipitation  intensity  and  change  and  to  note  any  consistent 

5.  "Wash,  C.  H.,  J.  Stremikis  and  D.  R.  Johnson,  1979:  Objective 
forecasts  of  subsynoptic  convective  areas  using  ■interactive 
computer  systems.  AMS  Preprint  Volume,  11th  Conference  on  Severe 
Local  Storms,  2-5  October  1979. 


tendencies  in  the  behavior  of  the  parameters. 

Only  precipitation  In  the  form  of  snow  was  considered  since  any 
intensity  changes  are  easily  detectable  from  changes  in  horizontal 
surface  visibility.  One  obvious  problem  is  the  influence  of  blowing 
snow  on  the  surface  visibility.  This  effect,  however,  was  considered 
less  critical  than  the  usual  presence  of  fog  and/or  other  obscurations 
in  continuous  rainfall  cases  which  makes  visibility  less  useful  as  an 
intensity  indicator.  Stations  which  exhibited  considerable  variability 
in  snowfall  intensity  (visibility)  were  selected  from  four  recent 
(1980-1981)  snowstorm  occurrences.  A  total  of  over  300  brightness 
arrays  centered  on  these  stations  formed  the  data  base  for  the  study. 
The  results  of  this  research,  with  addition  of  several  more  cases  and 
further  refinement  of  the  numerical  parameters  used,  could  aid  in 
development  of  techniques  for  detection  and  very  short  range 
forecasting  of  rapid  changes  in  snowfall  intensity  (visibility)  using 
only  GOES  data  and  surface  observations  for  information  sources. 


2.  BACKGROUND 


Since  the  advent  of  the  meteorological  satellite  in  the  late 
1950's,  research  has  been  directed  towards  the  pioneering  of  new 
methods  using  satellite  data  and  the  Improvement  of  older  methods. 

This  trend  continues  into  the  decade  of  the  1980's  with  large  research 
efforts  in  the  atmospheric  sciences  revolving  around  the  use  of 
satellite  Information-  With  the  launching  of  the  first  geosynchronous 
operational  satellite,  SMS- II,  In  May  of  1974,  a  breakthrough  in 
satellite  meteorology  occurred  with  the  capability  to  resolve  the  time 
evolution  of  atmospheric  phenomena.  Considerable  applications  have 
been  developed  for  the  use  of  GOES  data.  In  addition  to  diagnosis  and 
forecasting  through  viewing  imagery,  development  of  techniques  has  been 
completed  for: 

1)  determining  upper  level  wind  structure  from  observed  cloud 
motions6 

2)  estimating  convective  rainfall  rates,  from  IR  brightness  data,  in 
tropical  convection7 

3)  indirectly  sounding  the  temperature  and  moisture  structure  of  the 

'6~.  Johnson,  G.  L.,  and  D.  Suchman,  1980:  Intercomparisons  of  SMS  wind 
sets:  a  study  using  rapid  scan  imagery.  Mon.  Wea.  Rev.,  108, 
1672-1688. 

7.  Stout,  J.,  D.  Martin  and  D.  N.  Sikdar,  1979:  Estimating  rain  with 
geosynchronous  satellite  Images.  Mon.  Wea.  Rev.,  107,  585-598. 
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atmosphere8. 

The  work  conducted  in  this  study  is  most  closely  related  to  (2) 
except  that  the  research  deals  with  convection  in  mid-latitudes.  A 
somewhat  similar  study  conducted  by  Wylie9  attempted  to  apply 
techniques  developed  for  estimating  tropical  convective  precipitation 
rates10  to  middle  latitudes.  These  results  indicated  that,  although 
estimated  errors  for  these  techniques  were  much  larger  in  middle  than 
tropical  latitudes,  the  relative  accuracy  might  be  increased  with 
additional  parameterization  through  use  of  radiosonde  data.  Since  the 
techniques  were  originally  developed  for  convection  in  warm  air  masses, 
errors  in  the  application  of  the  tropical  techniques  to  mid-latitudes 
also  varied  greatly  between  the  summer  and  fall  cases  examined.  Little 
additional  work  has  been  conducted  either  to  apply  already  developed 
techniques  or  to  develop  new  techniques  for  determination  of 
mid-latitude  precipitation  solely  from  GOES  data. 


IT  Smith,  W.  L.,  V,  E.  Suomi,  W.  P.  Menzel ,  H.  M.  Woolf,  L.  A. 

Sromovsky,  H.  E.  Revercomb,  C.  M.  Hayden,  D.  N.  Erickson  and  F. 

R.  Mosher,  1981:  First  sounding  results  from  VAS-D.  Bull.  Amer. 
Meteor.  Soc.,  62,  232-236. 

9.  Griffeth,  C.  G.,  W.  L.  Woodley,  P.  G.  Grube,  D.  W.  Martin,  J.  Stout 

and  D.  N.  Sikdar,  1978:  Rain  estimation  from  geosynchronous 
satellite  imaqery-visible  and  infrared  studies.  Mon.  Wea.  Rev., 
106,  1153-1171. 

10.  Griffeth,  C.  G.,  W.  L.  Woodley,  P.  G.  Grube,  D.  W.  Martin,  J. 

Stout  and  D.  N.  Sikdar,  1978:  Rain  estimation  from 
geosynchronous  satellite  imagery-visible  and  infrared  studies. 
Mon.  Wea.  Rev. ,  106,  1153-1171. 
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Without  question,  a  large  amount  of  precipitation  in  mid-latitudes 
falls  as  continuous,  stratiform  type,  rain  or  snow.  Generally,  the 
intensity  of  this  type  of  precipitation  is  quite  uniform  with  little 
short-term  variation.  However,  convective  activity  is  occasionally 
imbedded  within  areas  of  stratiform  precipitation  and  under  these 
circumstances  very  large  temporal  changes  of  intensity  occur  with  the 
passage  of  convective  cells  over  a  particular  locale.  Often  these 
enhanced  precipitation  events  elude  the  synoptic  reporting  network  and, 
except  for  radar,  go  undetected.  The  capability  to  remotely  detect 
these  events  through  the  use  of  GOES  data  is  highly  desirable. 

Typically,  the  convection  imbedded  within  a  mid-latitude  storm 
system  is  located  in  the  cold  sector  of  the  cyclone  in  the  proximity  of 
the  surface  cyclone  center11.  It  appears  that  this  convective  activity 
is  primarily  associated  with  either  occluded  cyclones  or  rapidly 
developing  cyclones  along  the  eastern  coast  of  continents. 12  Figure  1 
shows  some  typical  storm  structures  and  the  convective  patterns 
frequently  observed  in  the  cold  sectors.  Note  the  banded  structure  of 
the  convection  commonly  associated  with  an  occlusion  and  the  cluster 
pattern  usually  observed  in  the  coastal  storms.  One  can  suggest  the 
following  reasons  for  the  preferred  convective  development  with  the 
occluded  and  coastal  cyclones. 

TH  Hobbs,  P.  V. ,  and  J.  D.  Locatelli,  1978:  Rainbands,  precipitation 
cores  a'd  generating  cells  in  a  cyclonic  storm.  J.  Atmos. 

Sci . ,  j5,  230-241. 

12.  Houze,  R.  A.,  Jr,,  P.  V.  Hobbs,  K.  R,  Biswas  and  W.  M.  Davis, 

1976:  Mesoscale  rainbands  in  extratropical  cyclones.  Mon.  Wea. 
Rev.,  104,  868-878. 


In  the  case  of  occluded  storms,  a  primary  forcing  for  the 
stratiform  precipitation  in  and  near  the  storm  center  is  the  ascent  due 
to  low  level  frictionally  induced  convergence.  With  the  occlusion,  a 
very  cold  pool  of  air  aloft  becomes  located  directly  over  the  surface 
cyclone  center  and  results  in  a  minimum  of  static  stability  in  the  mid 
troposphere  near  the  storm  center.  If  sufficient  frictional 
convergence  is  present  to  lift  the  moist  air  through  the  usually  stable 
lower  levels  of  the  atmosphere,  convective  activity  develops  in  the 
less  stable  layers  of  the  mid  troposphere.  Rapidly  deepening  coastal 
cyclones  are  nearly  always  associated  with  intense  Positive  V^orticity 
Advection  (PVA)  aloft  and  upward  vertical  motion.  This,  combined  with 
the  readily  available  moisture  supply  of  coastal  locations  and  thus 
relatively  high  equivalent  potential  temperatures,  likely  triggers 
convection  near  the  cyclone  center  where  upper  level  divergence  is 
maximized.  A  positive  feedback  results  with  the  strong  latent  heat 
release  forcing  the  enhancement  of  the  convection. 

The  convective  activity  just  described  usually  propagates  with  the 
middle  level  wind  flow  and  persists  for  a  period  of  several  hours.  As 
noted  by  previous  Investigators,  these  convective  areas  are  readily 
evident  in  visible  satellite  imagery  and  usually  are  associated  with 
bright  areas  in  the  IR  imagery.  Since  the  brightness  of  a  radiating 
body  in  an  IR  satellite  image,  as  conventionally  displayed,  is 
Inversely  proportional  to  the  absolute  temperature,  a  large  brightness 


value  would  Imply  a  cold  temperature13.  The  relative  high  brightness 
displayed  by  these  Imbedded  convective  areas  requires  that  the  cloud 
tops  of  the  convection  are  significantly  colder  than  the  tops  of  the 
surrounding  stratiform  cloudiness,  observations  from  aircraft  indicate 
that  the  cloud  tops  are  usually  several  thousands  of  feet  higher  than 
the  stratiform  cloud  tops14.  Due  to  the  fact  that  many  of  these 
convective  bands  and  clusters  are  very  small,  with  typical  radial 
dimensions  of  several  miles,  high  resolution  satellite  data  is  required 
to  detect  the  majority  of  the  convective  activity. 
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T3.  Brandi i,  H.  W. ,  1976:  Satellite  Meteorology,  Air  Weather  Service 
Tech.  Rep.,  AWS-Tl. -76-264,  Hq  Air  Weather  Service,  Scott  A^B, 

M  ,  f8S  pp. 

14.  Brandli,  op.  cit. 


3.  SELECTION  AND  DESCRIPTION  OF  SYNOPTIC  CASES 

For  the  screening  and  selection  of  the  SMS  satellite  information 
and  hourly  surface  observations  which  formed  the  d;.ta  base  for  this 
study,  fifteen  snowstorms  during  the  winters  of  1979-1980  and  1980-1981 
were  examined.  In  the  screening,  if  surface  visibility  observations 
for  several  stations  within  a  particular  storm  indicated  variable 
snowfall  intensity,  SMS-II  satellite  information  and  surface 
observations  were  selected  from  the  McIDAS  satellite  data  archive.  The 
imagery  with  2nm  resolution  at  subpoint  was  obtained  at  half-hour 
intervals  over  an  eight-hour  period  for  seven  of  the  storms  with 
variable  snowfall  intensity.  The  methods  of  obtaining  the 
station-centered  arrays  of  digital  IR  brightness  values  used  for  study 
and  a  detailed  description  of  their  structure  are  discussed  in  a  later 
section. 

Of  the  snowstorms  examined,  only  four  cases  satisfied  the 
requirement  that  the  actual  surface  observations  of  variable  visibility 
and  the  area  of  embedded  convection  evident  in  the  satellite  imagery 
were  co-located.  A  brief  synoptic  description  of  these  four  cases 
selected  for  detailed  study  is  now  presented. 

3.1  14  April  1980,  Central  United  States. 

The  Midwest  storm  on  14  April  1980  was  in  a  more  developed  stage 
than  any  of  the  other  three  cases  which  were  studied.  The  system  had 
occluded  prior  to  this  data  period. 
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The  Midwest  cyclone  of  14  April  1980  developed  along  the  Louisiana 
Gulf  coast  early  on  the  preceding  day.  At  that  time,  a  strong  low 
level  baroclinic  zone  was  situated  along  the  coast  and  a  very  high 
amplitude  trough  at  upper  levels  was  located  over  the  central  II.  S. 
with  an  axis  from  Lake  Winnepeg  in  southern  Canada  to  southern  Texas. 
Strong  westerlies  with  a  110  knot  wind  maximum  were  located  in  the 
southern  portion  of  the  trough.  The  patterns  of  divergence  and 
vorticity  indicated  that  the  atmospheric  structure  was  ideal  for  the 
occurrence  of  surface  cyclogenesis  in  the  central  Gulf  coast  states. 

The  surface  cyclone  deepened  during  its  northerly  movement  with 
the  central  sea  level  pressure  falling  8  mb  in  the  12  hour  period 
ending  00C0  GMT  14  April.  At  this  time,  the  storm  was  centered  in 
extreme  northern  Alabama.  Occlusion  was  occurring  rapidly  during  this 
period.  Based  on  analyses  from  the  National  Meteorological  Center 
(NMC),  the  storm  was  fully  occluded  by  1200  GMT  14  April  (Figure  2). 
During  the  following  IZ  hours,  the  center  continued  to  move  slowly 
northward  to  a  location  in  east  central  Illinois  by  0000  GMT  15  April. 
The  500  mb  flow  during  this  stage  of  occlusion,  1200  GMT  14  April  to 
0000  GMT  15  April,  was  dominated  by  an  expansive  cold  core  closed  low 
over  the  central  U.S.  The  circular  symmetry  of  the  height  contours 
shown  in  Figure  3  suggests  that  during  this  period  the  PVA  was  a 
minimum;  hence,  little  further  intensification  of  the  storm  occurred. 

The  significant  snowfall  during  14  April  was  confined  to  a  narrow 
band  from  northeastern  Arkansas  to  northern  lower  Michigan.  Figure  2a 
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indicates  the  stations  in  this  band  over  which  brightness  arrays  for 
the  study  were  centered.  The  embedded  convective  nature  of  the  clouds 
In  northern  Illinois  and  southern  Wisconsin  within  the  snow  band  is 
evident  from  the  satellite  Images  shown  in  Figures  5  and  6.  The 
convection  is  especially  apparent  in  the  visible  images.  Note  that  the 
large  brightness  values  in  the  IK  images  are  also  coincident  with  these 
areas.  Many  of  the  banded  convective  areas  persisted  for  several  hours 
during  their  relative  propagation  around  the  circulation  center. 

Strong  low  level  convergence  near  the  surface  cyclone  center  combined 
with  the  unstable  lapse  rate  shown  in  Figure  4  was  likely  a  primary 
factor  in  the  forcing  of  this  convective  activity. 

3.2  6  February  1980,  Southeastern  United  States. 

The  storm  in  the  Southeastern  United  States  on  6  February  1980  was 
a  type  very  different  rrom  the  Midwest  storm  that  was  just  described. 

At  the  beginning  of  the  satellite  data  period,  1400  GMT  6  February,  a 
weak  occluded  cyclone  that  had  been  drifting  slowly  southeastward  for 
several  days  accompanied  by  light  precipitation  was  located  in  the 
lower  Ohio  River  valley.  At  the  same  time,  with  favorable  upper-level 
and  lower-level  baroclinic  support  along  the  Atlantic  coast, 
cyclogenesis  was  occurring  along  the  South  Carolina  coast.  The  primary 
precipitation  area  was  associated  with  this  newly-developed  cyclone 
center. 

In  Figure  7,  the  development  of  the  surface  center  off  the  coast 
was  quite  slow  initially,  but  rapid  deepening  eventually  occurred  late 
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in  the  period  of  interest.  The  central  sea  level  pressure  dropped 
14  mb  in  the  six  hours  ending  0000  GMT  7  February.  Eventually,  the 
lowest  pressure  attained  by  the  storm  was  986  rub  at  1200  GMT 
7  February.  The  most  significant  snowfall  associated  with  the  new 
center  was  in  a  band  from  South  Carolina  northeastward  through 
Washington,  O.C. 

The  distribution  of  stations  reporting  variable  snowfall 
intensities  on  the  eastern  slopes  of  the  Appalachian  chain, 

(Figure  7a),  together  with  the  predominantly  easterly  low  level  wind 
flow,  suggests  that  orographical ly  induced  rising  motion  was  partially 
responsible  for  the  convective  activity  in  the  ccld  sector  of  this 
cyclone.  Unlike  the  first  case,  this  developing  system  was  linked  with 
a  ratner  strong,  fast-moving  upper-level  short  wave  disturbance.  The 
strong  PVA  over  the  region  of  development  indicated  in  the  500  mb  chart 
(Figure  8)  and  the  other  NMC  model  analyses  (not  shown)  likely  also 
aided  the  convective  development. 

The  visible  and  Irt  satellite  images  of  Figures  10  and  11  show 
considered1?  cirrus  cloudiness  present  throughout  the  data  period 
11400  "PT  -  2100  GMT  6  February).  The  convective  activity,  though 
somewhat  difficult  to  observe  in  visible  and  IR  channels  due  to  the 
presence  of  this  cirrus,  seemed  in  this  case  to  be  organized  in 
clusters  rather  than  the  banding  that  is  more  typical  of  occluded 
systems,  such  as  in  the  14  April  1980  case.  The  clusters  persisted  for 
the  entire  eight  hours  of  satellite  data.  Because  of  the  obscuring 
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cirrus  cloudiness,  the  detection  of  the  convective  clusters  with  the 
digital  IR  data  is  more  difficult  In  this  case  than  in  any  of  the 
others. 

3.3  2  March  1980,  Southeastern  United  States. 

The  storm  over  the  Southeastern  United  States  on  2  March  1980  Is  a 
classic  example  of  Cape  Hatteras  development.  At  0000  GMT  2  March,  a 
stationary  surface  frontal  system  was  located  just  off  the  coast  of  the 
U.  S.  Several  waves  of  low  pressure  were  situated  along  this  front 
with  the  dominant  center  being  located  approximately  400nm 
south-soutneast  of  Cape  Hatteras.  By  1200  GMT  2  March,  only  one  center 
situated  approximately  250nm  southeast  of  the  Cape  remained.  From  this 
time  through  the  end  of  the  data  period,  0000  GMT  3  March,  the  center 
intensified  rapidly  with  a  16  mb  central  pressure  drop  occurring  in  the 
twelve  hours  ending  0000  GMT  3  March.  A  comparison  of  Figures  12b  and 
12c  shows  how  slowly  the  cyclone  center  moved  during  the  period  of 
deepening. 

A  slight  northwest-southeast  orientation  of  the  500  mb  short  wave 
trough  axis,  acquired  by  0000  GMT  3  March,  suggests  very  intense  PVA 
just  east  of  the  wave  (Figure  13).  The  satellite  imagery  loops  suggest 
that  the  PVA  helped  initiate  the  convective  activity  associated  with 
this  storm.  As  the  PVA  encountered  the  coastal  frontal  zone  and  a 
rather  unstable  alrmass  overlying  the  frontal  inversion,  rapid 
development  of  small  convective  clusters  occurred.  Near  the 
termination  of  the  satellite  data  period,  2100  GMT  2  March,  an 
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extensive  cirrus  overcast  developed  over  and  to  the  east  of  the 
convective  activity  (Figure  16).  The  shield  did  not  interfere  with 
detection  of  the  convection  as  in  the  previous  case  because  of  Its 
advectlon  downstream.  Individual  convective  clusters  persisted  for 
several  hours  during  the  cyclone's  propagation  northeastward.  The 
convective  clusters  appealed  more  vigorous  in  the  imagery  loops  for 
this  case  than  any  of  the  other  case  studies. 

3.4  31  January  -  1  February  1981,  Central  United  States. 

The  fourth  and  final  case  In  the  study  was  a  rather  uncommon  one. 
The  heaviest  snowfall  area  was  not  associated  with  any  organized  storm 
center  but  rather  with  a  poorly  defined  north-south  oriented  surface 
trough  of  low  pressure.  Development  of  a  surface  cyclone  over  southern 
Kansas  was  just  beginning  by  0000  GMT  1  February  (Figure  17).  A  large 
anticyclone  of  Canadian  origin  had  moved  across  the  Great  Lakes  region 
during  30  January  and  was  centered  over  Ohio  by  1200  GMT  31  January. 
Associated  with  this  event,  troughing  at  the  surface  and  aloft  was 
occurring  over  the  Rocky  Mountains  and  western  Great  Plains.  The 
southerly  flow  between  the  trough  over  the  western  Great  Plains  and  the 
anticyclone  over  Ohio  was  quite  intense  and  extended  upward  to  at  least 
300  mb.  Extensive  warm  advection  was  occurring  through  a  deep  layer  of 
the  atmosphere  with  a  maximum  located  in  the  Nebraska-Kansas  region. 

The  broad  structure  and  northeast  to  southwest  orientation  of  the 
500  mb  trough  shown  in  Figure  18  does  not  suggest  strong  PVA  over  the 
central  Plains.  Hence,  warm  advectlon  was  apparently  the  primary 
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forcing  for  the  upward  motion  and  precipitation.  The  convective 
activity  was  embedded  within  the  region  of  maximum  warm  advectlon.  The 
convection  Increased  In  area  and  Intensity  during  the  afternoon  hours 
apparently  due  to  further  destabilization  above  the  750  mb  level 
through  differential  temperature  advectlon  (Figure  19).  However,  In 
this  case  all  the  bands  were  moving  toward  the  northeast  Instead  of 
spiraling  around  any  circulation  center.  Cirrus  cloudiness  was  present 
east  of  the  con\ect1ve  activity;  however,  none  was  noted  over  the 
region  of  convection  itself. 

During  the  day,  the  system  propagated  eastward  with  little 
Intensification.  After  0000  GMT  1  February,  the  system  intensified, 
accelerated  toward  the  northeast,  and  attained  a  central  surface 
pressure  of  993  mb  by  0000  GMT  2  February  over  the  James  Bay  region  of 
Canada. 


4.  DATA  SAMPLE 


The  hourly  surface  synoptic  reports  and  the  digital  1R  brightness 
arrays  used  In  this  research  were  obtained  through  Me  IDAS.  Data  for 
the  14  April,  6  February  and  2  March  1980  cases  were  retrieved  from  the 
system's  extensive  data  archive  containing  the  surface,  upper  air  and 
satellite  Information  received  since  approximately  1975.  The  data 
acquisition  for  the  final  case  was  completed  In  real  time.  In  each 
Instance,  &  listing  of  the  surface  observations  for  all  stations 
reporting  snowfall  of  any  Intensity  was  obtained  each  hour  for  a  period 
of  six  to  twelve  hours  (Figure  20a).  Next,  an  hour  by  hour  listing  was 
produced  for  any  station  which  reported  moderate  or  heavy  snowfall  (S 
or  S+)  for  at  least  one  hour.  An  example  for  Ralelgh-Durham,  North 
Carolina  is  shown  In  Figure  20b. 

The  selection  of  location  and  time  for  which  to  obtain 
corresponding  IR  brightness  arrays  was  based  on  the  hourly  variation  of 
visibility  for  this  subset  of  stations.  A  square  cursor,  defined  to  be 
18  television  pixels  on  a  side,  was  centered  on  each  of  the  desired 
stations  in  the  hourly  image.  An  array  of  visible  and  IR  brightness 
values  ranging  between  0  and  255  were  produced,  an  example  of  which  Is 
given  In  Figure  21a.  Each  array  contains  324  data  points,  and  each 
data  point  represents  a  2nm  x  2nm  areal  average  brightness  value.  The 
total  horizontal  area  encompassed  by  each  array  is  approximately  1300 
square  nautical  miles  with  small  differences  depending  on  the  latitude 
and  longitude  of  the  centrally  located  station.  From  each  array  the 
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mean,  median,  mode  and  standard  deviation  of  the  brightness  values  were 
also  obtained  (Figure  21b).  A  complete  listing  of  stations  and  hours 
for  which  arrays  were  produced  in  each  case  Is  shown  in  Tables  1 
through  4. 
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tout  (on) 


STATION' 

it 

17 

IB 

It 

20 

21 

22 

21 

Kidtion,  WI 

z 

X 

X 

X 

Creen  lay^  UI 

X 

X 

X 

Sprlnffleld,  IL 

X 

X 

X 

Hanlcovoc,  VI 

X 

X 

Norli.  IL 

Z 

X 

OlCienva.  IA 

z 

X 

X 

X 

X 

X 

KoUni.  XL 

X 

X 

X 

X 

X 

X 

Dubuque,  XA 

z 

X 

X 

X 

Al p«n  * ,  HI 

X 

X 

X 

X 

Table  1:  Stations  and  hours  for  which  data  arrays  were  produced  for  14 
April  1980  case. 


HOUt  (Oft) 


STATION 

it 

IS 

It 

27 

la 

It 

20 

21 

lUlelith-Durhaa.  NC 

X 

X 

X 

X 

X 

Charlotte.  NC 

X 

X 

X 

X 

Columbia,  SC 

X 

X 

LrnchburR.  VA 

X 

X 

X 

X 

X 

X 

X 

Roanoke,  VA 

X 

X 

X 

X 

X 

Danville.  VA 

X 

X 

X 

Creaniboro,  NC 

X 

X 

Dull#*  Xnt.  Arpc. 

X 

X 

X 

X 

X 

Washington  Nat.  Arpt. 

X 

X 

X 

HartlnabuiR,  W 

z 

X 

X 

X 

X 

X 

Baltlrore,  KD 

X 

X 

X 

X 

Ulnston-Salta,  NC 

X 

z 

X 

CraanvUle.  SC 

X 

X 

lewlaburt,  VV 

X 

X 

X 

X 

X 

X 

X 

X 

Table  2:  Stations  and  hours  for  which  data  arrays  were  produced  for  6 
February  1980  case. 
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HOURS  (CKT) 


Hyannla  ,  HA _ 

Seymour,  HC _ 

Danville.  VA _ X_ 

Illiabeth  CltT  ,  NC  X 

Pope  A/ b .  WC _ 

Octane  HAS,  VA _ 

Hot  folk  .  VA _ t_ 

Langley  AT 8 ,  VA _ X 

AnjtrioWj  SC _ 

RockyHount .  SC _ X 

Hickory .  HC _ X_ 

Albany.  CA  X 


Fayet  t evil le .  NC 


Ralt lgh-Durhea ,  NC  X 


Atlanta.  CA _ 

Wlnaton-S a l an.  NC 


Cratnaboro,  HC 
Nantucket .  HA 


Table  3:  Stations  and  hours  for  which  data  arrays  were  produced  for  2 
March  1980  case. 


HOUl  {(XT) 


Norfolk,  HE 
Sidney,  HC 


(Lana a ■  Ctty  Int  .  Arpt. 


VLanaaa  City  Muni.  Arpt . 


Hlbblng,  W 
Pea  Holpca,  U_ 


Topeka,  C 


St.  Joaephe.  HO 
Mullen.  HE 


Kanhattan,  KS 
Souls  Cite.  1A 


Beatrice,  ME 
Lincoln.  WE 


McCook.  HE 


Alliance.  HE 
Alnavorth,  HC 


Table  4:  Stations  and  hours  for  which  data  arrays  were  produced  for  31 
January-1  February  1981  case. 


5.  PROCEDURES  FOR  CALCULATION  OF  PARAMETERS 


In  this  study  of  the  inference  of  convective  precipitation 
intensity  and  intensity  changes  from  geosynchronous  satellite 
information,  the  report  examines  the  following  parameters: 

1)  deviation  of  the  station  point  brightness  value  from  the  array 
areal  average  brightness,  and 

2)  the  Laplacian  of  the  brightness  field  estimated  by  four 
different  methods. 

The  following  procedures  were  applied  to  each  of  the  four  cases 
investigated: 

5.1  Selection  of  appropriate  brightness  arrays. 

For  a  brightness  array  to  be  selected,  (1)  current  hour  snowfall 
intensity  at  the  station  about  which  the  array  is-  centered  must  differ 
from  the  intensity  for  either  the  previous  or  following  hour  and  (2) 
standard  deviation  of  the  brightness  values  in  the  array  must  equal  or 
exceed  2.0.  Hence,  only  instances  which  displayed  hourly  intensity 
variations  as  well  as  variable  cloud  structure  in  the  IR  data  were 
employed  in  the  study.  These  restrictions  reduced  the  number  of  arrays 
used  for  computations  to  approximately  200,  from  more  than  300  arrays 
originally  obtained  from  the  four  cases.  A  subset  of  50  of  these 
arrays  which  exhibited  the  most  extreme  variations  in  snowfall  and 
cloud  brightnesses  (standard  deviation  >6.0)  was  considered  separately, 
and  a  comparison  of  results  is  made  between  this  subset  and  the  larger 
array  data  base  in  the  results. 
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5.2  Assignment  of  point  brightness  values. 

Although  McIDAS  misalignment  In  point  navigation  routines  is 
assumed  to  be  no  greater  than  lnm  for  SMS-II  data)b  all  brightness 
values  used  in  calculations  were  obtained  by  averaging  surrounding 
values  to  reduce  navigation  biases.  The  average  poin.  brightness  was 
calculated  from  four  surrounding  points  located  equidistant  from  the 
central  point.  An  example  of  the  average  point  brightness  calculation 
is  presented  in  Figure  22,  where  the  point  brightness  for  the  station 
at  X  is  obtained  by  averaging  the  four  surrounding  values.  The  values 
on  the  array  perimeter  were  not  averaged. 

5.3  Calculation  of  the  five  parameters  for  selected  arrays. 

In  this  study,  the  first  parameter  computed  from  each  array  was 
the  deviation  of  the  station  brightness  (DM)  from  the  mean  areal 
average  brightness  of  the  array.  Positive  values  for  DM  imply  that  the 
cloud  top  above  the  station  is  colder  and  higher  than  the  average  cloud 
top  in  the  array  field.  Similarly,  negative  values  indicate  a  warmer 
and  lower  cloud  top  than  the  average  value. 

The  two-dimensional  Laplacian  values  of  the  brightness  field  cen¬ 
tered  on  the  station  were  the  second  parameters  computed.  The  two- 
dimensional  Laplacian  operator  v2 

TIT!  Griffeth,  C.  G.,  W.  L.  Woodley,  P.  G.  Grube,  D.  W.  Martin,  J. 

Stout  and  D.  N.  Sikdar,  1978:  Rain  estimation  from 
geosynchronous  satellite  imagery-visible  and  infrared  studies. 
Mon.  Wea.  Rev.,  106,  1153-1171. 
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32  32 

V2  *  —  +  --- 
3x2  3y2 

in  finite  difference  form  is 

2  (  )l  +  (  )2  +  (  )3  +  {  )4  "  4(  )0 

7°" . ix2 . 

Here  the  subscripts  indicate  evaluation  from  scalar  values,  at  the 
discrete  points  of  the  sketched  grid  array, 

2 

1 - o - 3 

i 

I 

4 

In  this  configuration  to  estimate  the  Laplacian,  four  of  the  points  are 
distributed  equidistant  and  orthogonal  around  a  central  point.  The 
two-dimensional  Laplacian  of  any  scalar  field  is  simply  a  measure  of 
the  rate  of  change  of  the  horizontal  gradient  on  a  horizontal  plane, 
and  is  used  to  locate  both  relative  maxima  and  minima  of  values  in  the 
scalar  field.  Its  sign  is  positive  for  relative  minima  and  negative 
for  relative  maxima.  When  applied  to  an  IR  brightness  field,  the 
Laplacian  effectively  indicates  the  increased  (decreased)  vertical 
extent  of  a  cloud  top  over  the  station,  relative  to  the  cloud  tops  at 
the  perimeter  points,  associated  with  a  brightness  maxima  (minima). 

This  is  accomplished  through  the  cloud  top  height-temperature-IR 
brightness  relationship  previously  discussed.  As  the  "cl rcularlty"  of 
these  localized  maxima  and  minima  of  brightness  increase,  so  does  the 
magnitude  of  the  Laplacian.  An  example  of  the  dependence  of  magnitude 
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upon  "circularity"  Is  shown  in  Figure  23. 

The  computed  five  point  Laplacian  (LAP  1 )  values  were  highly 
dependent  upon  the  orientation  of  the  bright  IR  bands  with  respect  to 
the  orthogonal  axes  along  which  the  discrete  points  are  located. 

Figure  24  illustrates  a  situation  where  the  five-point  method  results 
are  misleading,  due  to  dependence  upon  orientation.  This  dependency  is 
a  problem  for  the  convective  activity  in  the  banded  structure  of  cases 
one  and  four. 

In  an  attempt  to  limit  this  problem  a  Laplacian  with  a  nine-point 
structure  (LAP2)  based  on  a  polar  coordinate  geometry  was  selected. 
Eight  points  are  located  equidistantly  and  symmetrically  around  a 
central  point.  Both  the  nine-  and  five-point  structures  behave 
similarly  In  terms  of  sign  and  magnitude  characteristics.  See  Figure 
25  for  schematic  representations  of  both  Laplacian  structures  along 
with  the  equations  used  for  calculations.  The  subscripted  B  values 
represent  brightnesses  at  the  respective  numbered  points  used  in  the 
calculations. 

The  LAP1  and  LAP2  were  each  calculated  twice  for  all  arrays;  the 
ones  with  a  grid  interval  of  approximately  8  nm  to  be  called  the  inner 
Laplacians  (ILAP1  and  ILAP2),  and  the  others  with  an  interval  of  16  nm 
called  the  outer  Laplacians  (0LAP1  and  0LAP2).  Intercomparisons 
between  Laplacians  with  different  grid  intervals  were  not  made,  and 
only  the  sign  and  relative  magnitude  of  each  type  was  considered. 

Hence,  division  by  the  square  of  the  point  separation  (aX2)  was  not 
necessary. 
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5.4  Parameter  categorization. 

The  resultant  values  of  the  computation  of  the  five  parameters 
(DM,  ILAP1,  ILAP2,  0LAP1,  0LAP2),  were  compiled  for  individual  stations 
and  grouped  by  case  number  in  tables  similar  to  the  example  in  Figure 
26  for  Madison,  Wisconsin  for  the  14  April  1980  data,  case  number  1. 

The  visibility  (VIS)  and  weather  (WX)  data  were  obtained  from  the  sta¬ 
tion  listings  described  earlier.  From  the  tabulated  data,  frequency 
distri butions  and  mean  values  for  the  five  parameters  were  generated 
and  categorized  according  to  snowfall  intensity.  Two  categories  were 
considered,  light  snow  (VIS>0.7  miles)  and  moderate  or  heavy  snowfall 
( VIS<0. 7  miles).  Potential  relationships  between  the  five  parameters 
and  the  intensity  of  snowfall  as  it  varied  over  space,  from  station  to 
station,  are  scrutinized  in  the  following  section. 

Another  important  consideration  of  this  study  was  the  temporal  be¬ 
havior  of  each  parameter  at  a  single  location.  Changes  of  each  of  the 
five  parameters  during  one  hour  periods  were  calculated  for  all  instan¬ 
ces  where  snowfall  intensity  changed  during  the  preceding  hour  (ADM, 
AILAP1,  AILAP2,  aOLAPI,  A0LAP2).  Here  the  resultant  values  are  cate¬ 
gorized  under  one  of  these  headings: 

1)  positive  intensity  change  (visibility  reduction), 

2)  negative  intensity  change  (visibility  increase). 

Hourly  changes  of  the  five  parameters,  as  well  as  visibility,  computed 
from  the  data  in  Figure  26  were  compiled  for  each  synoptic  case  in 
tables  headed  according  to  intensity  change.  Figure  27  shows  an  exam¬ 
ple  of  the  tabulation  format  for  intensity  increases. 
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6.  RESULTS 

The  results  for  estimates  of  precipitation  intensity  and  Intensity 
changes  from  GOES  digital  IR  brightness  for  four  rises  of  middle 
latitude  convective  activity  imbedded  in  continuous  stratiform 
precipitation  are  now  presented.  First,  a  comparison  is  made  between 
means  and  frequency  histograms  of  each  parameter  for  the  categories  of 
spatially  and  temporally  varying  snowfall  intensity,  (S-)  and  (S,S+). 
The  hourly  changes  in  each  parameter,  means  and  frequency  plots,  are 
then  computed  and  their  ability  to  determine  cases  of  intensity 
increases  and  decreases  with  time  is  tested.  Next,  an  analysis  of  the 
most  skillful  indicators  of  snowfall  intensity  and  intensity  change  is 
made.  Finally,  the  above  comparisons  for  intensity  change  are  made  for 
the  special  50  array  subset  with  standard  deviation  of  brightness 
values  equal  to  or  exceeding  6.0. 

6.1  Comparison  of  snowfall  parameters  with  observed  snowfall  intensity. 

Of  the  210  IR  brightness  arrays  meeting  the  mandatory  requirements 
set  forth  earlier,  136  were  centered  on  stations  reporting  light 
snowfall,  while  the  remaining  74  arrays  were  for  stations  reporting 
moderate  or  heavy  snowfall.  Average  values  of  all  five  parameters  and 
visibility  were  calculated  for  these  two  intensity  categories,  and  the 
results  are  presented  in  Table  5.  All  of  the  tests  are  based  on  the 
hypothesis  that,  for  convection  embedded  in  regions  of  continuous  snow, 
a  station  reporting  a  heavy  snowfall  intensity  would  have  a  cloud  top 
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Table  5:  Average  values  of  all  parameters  and  visibility  for  each 
snowfall  intensity  category. 
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of  much  larger  vertical  extent  located  above  It,  relative  to  nearby 
areas.  Thus,  one  would  expect  a  relative  maximum  of  IR  brightness  In 
the  array  at  the  station  location.  The  distribution  of  the  signs  of 
the  computed  parameters  shown  In  Table  6  idealized  by  this  hypothesis 
is  likely  valid  for  the  cases  examined  when  the  areal  extent  of 
relative  maxima  or  minima  of  intensity  and  brightness  are  encompassed 
within  the  36  nm  x  36  nm  area  of  the  array. 

Mean  values  of  all  the  parameters  display  the  expected  sign  for 
the  light  snow  category.  However,  the  sign  expectations  for  the 
heavier  (S,S+)  occurrences  were  only  satisfied  for  0LAP1,  and  even  then 
its  value  is  only  slightly  negative.  The  values  in  Table  5  suggest 
..hat  an  occurrence  of  light  snowfall  with  an  average  visibility  of  1.7 
miles  is  indeed  associated  with  an  IR  brightness  minimum  relative  to 
nearby  areas,  but  a  similar  interpretation  cannot  be  made  for  the 
heavier  intensity  occurrences  from  the  four  cases. 

Mean  parameter  values  for  each  of  the  synoptic  cases  individually 
are  shown  in  Table  7.  Recall  that  in  cases  one  and  four  the  imbedded 
convection  displayed  a  banded  structure,  while  in  cases  two  and  three 
the  convection  is  cellular.  The  computational  results  in  Table  7  for 
the  nine  point  Laplacian  averages,  both  small  and  large  scale,  for 
cases  one  and  four  resulted  in  the  hyopthesized  sign  distribution. 

Such  was  not  the  situation  for  cases  two  and  three  involving  primarily 
cellular  convection  (Table  7). 

The  frequency  histograms  for  all  five  parameters  in  Figures  28 
through  30  exhibit  large  variance  of  values  for  both  intensity  cate- 
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Table  7:  Parameter  averages  and  number  of  occurrences  for  each 
snowfall  intensity  category  for  individual  synoptic  cases. 
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gorles.  Standard  deviations  for  each  are  shown  in  Table  8,  The  large 
variances,  coupled  with  the  near  zero  values  for  all  parameter 
averages,  Imply  that  within  the  array  the  overall  relationship  between 
the  means  of  the  parameter  and  snowfall  intensity  is  difficult  to 
determine.  However,  a  relationship  Is  seen  between  the  nine-point 
Laplacian  averages  and  intensity  for  cases  of  banded  convection. 

One  subtle  difference  between  the  Intensity  categories  for  each 
parameter  is  in  the  skewness  of  the  frequency  distributions.  Positive 
skewness  implies  a  tendency  for  observations  to  cluster  In  the  "pper 
half  of  the  range  of  values,  while  negative  skewness  means  that  more 
observations  are  found  in  the  lower  half.  Figures  28  through  30 
present  frequency  histograms  of  the  DM,  ILAP1,  0LAP1,  ILAP2,  0LAP2 
values  for  both  Intensity  categories.  Note  that  distribution  of 
parameters  within  each  frequency  interval  is  given  by  fj  <  f  <  fj  +  Af. 

Little  skewness  is  evident  In  the  distribution  of  station 
brightness  deviation  (DM)  from  the  array  mean,  but  rather  a  flat 
Gaussian  type  distribution  with  a  maximum  frequency  near  zero  is  found. 

The  distributions  for  the  two  outer  Laplacians  (0LAP1  and  0LAP2) 
display  slight  positive  skewness  for  light  snowfall  occurrences,  and  a 
negative  skewness  for  the  heavier  intensities.  This  suggests  a 
relative  brightness  maximum  is  associated  with  the  majority  of  heavy 
(S,S+)  occurrences,  and  a  relative  minimum  for  the  light  intensities. 
However,  both  inner  Laplacian  distributions  have  a  maximum  frequency  in 
the  interval  immediately  above  zero  for  the  heavier  intensities 


suggesting  relative  minima  of  brightness  are  usually  associated  with 
the  heavy  snowfall  occurrences,  a  result  opposite  to  the  hypothesized 
distribution. 

In  summary,  the  0LAP1  and  0LAP2  distributions  are  skewed  somewhat 
in  the  expected  directions,  and  of  those  two,  the  nine-point  structure 
(0LAP2)  seems  to  produce  the  most  favorable  results  when  compared  to 
the  hypothesized  brightness-intensity  relationship  presented  earlier. 
However,  a  definite  Inference  of  snowfall  intensity,  based  solely  upon 
the  five  parameters  just  examined,  does  not  seem  promising. 

One  possible  reason  for  the  failure  of  these  methods  to  adequately 

detect  the  regions  of  enhanced  snowfall  intensity  is  a  horizontal 

displacement  between  the  heavy  snowfall  area  at  the  surface  and  the 

convective  clouds  producing  the  heavy  snow.  This  displacement  may 

occur  with  strong  wind  shear,  in  both  direction  and  speed,  which  cause 

the  snow  to  blow  a  considerable  horizontal  distance  from  its  point  of 

formation  in  the  clouds  until  it  reaches  the  earth's  surface.  With  a 

4  km  level  of  origin  and  a  1  m/s  fall  velocity16,  a  mean  wind 

of  4  m/sec  relative  to  the  cloud  of  origin  and  the  earth's  surface  will 

produce  a  typical  displacement  of  12-16  km.  The  mean  wind  in  this  case 

is  defined  to  be  the  average  wind  between  the  earth's  surface  and  the 

level  of  origin.  For  a  linear  increase  with  height,  this  relative  mean 

wind  of  4  m/sec  corresponds  to  a  wind  difference  of  8  m/sec  from  the 

earth's  surface  to  the  level  of  origin.  Strong  wind  shear  was  present 

Tff.  Rogers,  R.  R.,  1979:  A  Short  Course  in  Cloud  Physics.  Pergamon 
International  Library,  235  pp. 
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in  all  cases  examined.  This  effect  is  likely  enhanced  for  very  narrow 
bands  or  clusters  of  enhanced  intensity. 

Figure  31  shows  the  Madison,  Wisconsin  centered  IR  brightness  from 
Case  1  arrays  for  1900,  2000  and  2100  GMT,  contoured  every  10  units, 

superimposed  on  the  concurrent  WSR57  radar  depictions  from  Neenah, 

\ 

Wisconsin.  The  shaded  areas  represent  radar  returns  of  level  two  or 
greater  (scaled  from  1  to  6).  Note  the  alignment  of  the  heavier 
precipitation  areas,  inside  the  array  boundaries,  with  the  largest 
brightness  values.  At  1900  GMT  and  again  at  2000  GMT,  MSN  reported  0.5 
mile  visibility  in  moderate  snowfall,  even  though  the  orientation  of 
brightness  maxima  within  the  array  would  imply  a  relative  minimum  of 
brightness  directly  over  the  station  itself.  The  largest  values  of 
both  IR  brightness  and  radar  return  were  located  to  the  east-northeast 
of  MSN.  Wind  flow  at  all  levels  up  to  and  including  300  mb  was  from 
the  northeast  at  1200  GMT  with  the  maximum  speed  of  26  m/sec  at  850  mb. 
This  suggests  that  the  snow  did  indeed  blow  southwestward  as  it  fell, 
displacing  the  region  of  heavy  snowfall  at  the  surface  from  the  clouds 
which  produced  it. 

The  problem  of  downwind  displacement  of  heavy  snow  due  to  wind 
shear  was  not  anticipated  in  the  initiation  of  this  research.  However, 
during  the  conduct  of  the  research,  facets  of  this  problem  became 
evident.  The  more  detailed  analysis  including  radar  information  from 
case  1  suggests  that  displacement  of  the  snow  between  the  convective 
clouds  of  origin  and  the  point  of  deposition  is  a  likely  reason  for  the 
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lack  of  success  In  the  statistical  analysis.  It  seems  that  the 
displacement  problem  could  be  reduced  by  using  a  technique  which 
centers  the  arrays  of  long  lived,  easily  tracked  brightness  maxima  and 
minima.  A  comparison  of  parameter  values  for  those  stations  reporting 
variable  snowfall  intensities  near  and  downwind  of  the  array  center 
could  be  conducted.  A  comparison  of  results  from  such  a  study  with  the 
results  of  this  research  would  be  most  informative. 

6.2  Comparison  of  snowfall  parameter  changes  with  observed  hourly 

snowfall  intensity  changes. 

During  the  four  cases  a  total  of  55  intensity  increases  and  62 
decreases  occurred.  In  these  events  hourly  changes  for  each  of  the 
five  parameters  and  visibility  were  computed  from  the  data  and 
tabulated  as  described  earlier  (Figure  27).  The  mean  values  are  given 
in  Table  9. 

Just  as  a  relative  maximum  in  the  IR  brightness  field  over 
stations  reporting  moderate  to  heavy  snowfall  would  be  expected,  the 
hypothesis  anticipates  that  an  increase  in  relative  brightness  would 
ideally  be  found  over  stations  experiencing  an  increase  in  intensity. 
Likewise,  a  decrease  in  relative  brightness  values  would  be  found  over 
stations  reporting  decreases  in  snowfall  intensity.  If  this  were 
indeed  the  case,  the  sign  distribution  seen  in  Table  10  would  be 
expected  for  the  parameter  change  averages  for  each  intensity  change 
category.  This  expected  distribution  was,  in  fact,  realized  for  all 
but  two  of  the  ten  average  parameter  change-intensity  change 
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Table  9:  Average  hourly  changes  In  all  parameters  for  each  intensity 
change  category. 
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Table  10:  Hypothesized  sign  distribution  of  average  hourly  parameter 
changes  for  each  intensity  change  category. 
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combinations,  ADM  and  TTDTPT.  Magnitudes  of  the  average  changes, 
of  all  parameters  except  for  &0W,  were  larger  than  the  corresponding 
magnitudes  of  the  instantaneous  values  for  each  hour  seen  in  Table  5 
earlier. 

Figures  32  through  34  contain  the  frequency  histograms  of 
parameter  changes  for  both  intensity  change  categories.  The  aDM  and 
aILAPI  distributions  show  very  little,  if  any,  skewness,  with  maximum 
frequencies  very  near  zero.  Of  the  remaining  distributions,  the 
nine-point  Laplacian  changes  show  the  largest  negative  shift  in  their 
distributions  for  intensity  increases,  and  slight  positive  shifts  for 
decreases  in  intensity.  In  particular,  a  considerable  negative 
displacement  of  the  frequency  maximum  is  noted  for  intensity  increases 
in  the  distributions  of  AILAP2  and  A0LAP2  (Figure  3>t).  This  accurately 
represents  the  hypothesized  distributions,  and  suggests  that  snowfall 
intensity  increases  might  well  be  indicated  by  hourly  decreases  in 
values  of  the  nine  point  Laplacians.  In  general,  the  skewness 
displayed  by  the  parameter  change-intensity  change  distributions  is 
larger  than  that  for  the  parameter  value-intensity  histograms  shown 
earlier.  Also,  when  oresent,  the  skewness  is  in  the  expected 
d i recti  on. 

6.3  Oetermi nation  of  most  skillful  indicators  for  intensity  and 

intensity  change. 

The  ability  of  the  five  parameters  and  their  charges  to 
discriminate  between  each  Intensity  and  intensity  change  category  was 
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tested  through  computation  of  a  frequency  difference  for  each  Interval. 
This  was  accomplished  by  subtracting  the  (S,S+)  frequency  value  from 
the  corresponding  (S-)  frequency  value  In  a  given  interval  for  a 
particular  parameter.  The  same  procedures  were  applied  to  the  hourly 
change  histograms,  with  the  Intensity  increase  frequency  in  an  interval 
subtracted  from  the  corresponding  Intensity  decrease  frequency  for  each 
parameter  change. 

Based  upon  the  idealized  relationships  of  IR  brightness  to 
Intensity  and  brightness  change  to  Intensity  change  discussed  earlier, 
an  (S-)  bias  In  the  difference  of  the  frequencies  is  expected  for  DM 
values  less  than  zero  and  a  (S,S+)  bias  in  the  difference  for  positive 
DM  values.  Similarly,  the  frequency  difference  of  Laplacian  values 
should  favor  (S-)  in  positive  intervals  and  (S,S+)  in  negative 
Intervals.  The  frequency  difference  computations  are  a  measure  of 
histogram  overlap  between  the  two  intensity  and  Intensity  change 
categories  considered  for  each  parameter.  A  perfect  overlap  (two 
identical  distributions)  would  have  a  zero  difference  for  each 
frequency  interval,  would  be  indicative  of  random  distributions  for  the 
snowfall  parameter  and  would  imply  no  capability  to  discriminate 
snowfall  intensity  at  the  earth's  surface.  Conversely,  no  overlap  at 
all  (two  completely  dissimilar  distributions)  would  indicate  certainty 
in  discriminating  snowfall  intensity  or  change. 

The  frequency  difference  distributions  for  each  parameter  are 
shown  In  Figures  35  and  36.  In  these  figures,  the  subtotal  for 
positive  and  negative  intervals  of  each  parameter  and  parameter  change 
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are  plotted  in  the  appropriate  quadrant  of  each  diagram.  In  Figure  35 
for  example,  the  difference  frequencies  for  DM  show  a  subtotal  for 
positive  Intervals  of  one  percent  In  favor  of  (S-)  occurrences,  while 
the  subtotal  frequency  difference  for  the  negative  Intervals  favors 
(S,S+)  occurrences  by  5  per  cent.  These  subtotals  are  extremely  small 
for  both  intervals.  In  fact,  the  distributions  for  (S-)  and  (S,S+) 
observations  are  nearly  perfectly  overlapped;  hence,  inference  of 
intensity  using  DM  as  an  indicator  was  not  justified.  In  the  four 
cases  examined,  any  randomly  chosen  DM  value  for  the  array  centered 
station  would  have  a  nearly  equal  likelihood  of  being  associated  with 
an  (S-)  or  an  (S,S+)  observation. 

Only  the  frequency  differences  for  0LAP1  and  0LAP2  display  the 
expected  distribution,  with  the  largest  sum  of  differences  being  22  for 
0LAP1.  Based  upon  the  data  from  the  synoptic  cases  chosen,  these 
results  indicate  that  if  negative  values  of  0LAP1  were  obtained  from 
the  array  data  base,  only  10  per  cent  would  be  uniquely  associated  with 
(S,S+)  observations  while  90  per  cent  would  have  been  associated  with 
either  (S-)  or  (S,S+)  occurrences.  Similarly,  only  12  percent  of  the 
positive  values  obtained  would  be  uniquely  associated  with  (S-) 
occurrences  and  88  percent  could  be  associated  with  either  intensity. 
Thus,  an  accurate  determination  of  snowfall  intensity  would  have  been 
possible  ^or  only  11  percent  of  the  total  0LAP1  values  obtained  under 
the  assumption  of  an  equal  probability  for  positive  and  negative  values 
in  the  sample.  Since  OLAP  appears  to  be  the  most  skillful  intensity 
indicator  of  the  five  parameters  considered,  it  was  concluded  that 
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Inference  of  snowfall  intensity  using  the  methods  of  this  study  where 
the  centroid  of  the  Laplacian  snowfall  parameter  is  located  with  the 
observing  station  does  not  yield  reasonable  accuracy.  However,  the 
colocation  of  areas  of  radar  return  and  brightness  values  in  the  case  1 
suggests  that  discriminiation  would  be  sharpened  if  allowances  were 
made  for  the  displacement  of  falling  snow  between  level  of  origin  and 
surface  deposition  due  to  wind  shear. 

The  distributions  of  frequency  differences  for  all  parameter 
changes  are  presented  in  Figures  37  and  38.  AILAP2  and  A0LAP2  have  the 
largest  subtotal  as  well  as  total  frequency  differences.  Correct 
determination  of  intensity  change  was  possible  for  27  percent  and  29 
percent  of  the  AILAP2  and  A0LAP2  samples,  respectively.  This  is  a 
significant  percentage  of  intensity  changes.  Further  work  in 
development  and  modification  of  the  methods  used  in  this  study  in 
conjunction  with  snowfall  displacement  due  to  wind  shear  should  be 
undertaken  to  improve  the  reliability  of  intensity  change  inferences. 

6.4  Determining  intensity  change  for  the  50  array  subset  (standard 
deviation  of  brightness  in  the  array  of  at  least  6.0). 

The  selection  of  a  restricted  data  subset  with  a  standard 
deviation  of  at  least  5.0  for  array  values  means  that  only  those  arrays 
with  extreme  variations  of  IR  brightness  are  considered.  The 
magnitudes  of  the  parameter  changes  for  this  subset  (Table  11)  are 
consequently  larger,  than  are  the  average  magnitudes  for  parameter 
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„ble  11:  Average  hourly  parameter  changes  for  each  intensity  change 
category,  for  the  special  50  array  subset  (standard 
deviation  of  brightness  values  >6.0). 
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changes  for  the  entire  sample.  The  total  number  of  Intensity  changes 
In  each  category  with  31  increases  and  30  decreases  is  approximately 
half  that  for  the  entire  array  data  set. 

Within  this  subset,  all  parameter  change  averages  display  the 
expected  sign  (Table  10).  The  frequency  histograms  also  exhibit  more 
displacement  between  maximum  frequency  intervals  for  cases  of  intensity 
increase  versus  decrease  (Figures  39  through  41).  However,  the 
distribution  overlap  is  also  greater,  thus  the  reliability  of  intensity 
change  determination  inferred  from  the  frequency  difference  diagrams, 
of  Figures  42  and  43,  is  not  increased.  The  most  skillful  indir  of 
intensity  change  for  the  50  array  subset  intensity  change  is  AILAP2 
with  a  total  frequency  difference  of  57,  less  than  the  most  skillful 
indicator  of  intensity  change  for  the  entire  array  set. 

In  this  array  subset,  large  standard  deviations  of  brightness 
values  of  small  scale  horizontal  variations  In  cloud  structure  are 
likely  present.  Thus,  the  ability  of  the  Laplacian  to  accurately 
resolve  smaller  scale  maxima  and  minima  in  the  brightness  fields  may  be 
decreased.  In  an  extreme  case,  possibly  several  maxima  and/or  minima 
may  exist  within  the  finite  grid  area  used  for  the  Laplacian  estimates 
within  a  given  array.  Under  these  circumstances ,  the  discrete 
estimates  of  the  Laplacians  of  the  brightness  data  may  indicate 
erroneous  values.  To  analyze  this  problem,  higher  resolution  satellite 
imagery  than  that  produced  for  this  study  must  be  used.  Also,  the 
problem  of  displacement  of  the  snow  by  wind  shear  may  be  compounded  for 
very  small  scale  convection. 
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7.  CONCLUSIONS 


This  report  presents  a  pilot  study  designed  to  develop  techniques 
applicable  to  middle  latitudes  for  inferring  precipitation  intensity 
and  intensity  changes  in  cases  of  continuous  precipitation  from  GOES  IR 
brightness  data.  Although  the  results  obtained  in  this  study  involved 
precipitation  only  in  the  form  of  snow,  they  likely  also  are  applicable 
to  studies  of  intensity  and  intensity  change  of  showers  embedded  in 
continuous  rain.  Snow  was  chosen  because  of  the  close  relation  between 
intensity  and  horizontal  surface  visibility. 

The  satellite  and  conventional  data  for  use  in  the  study  ware 
obtained  through  McIDAS  at  the  University  of  Wisconsin-Madison.  A  data 
base  consisting  of  hourly  digital  IR  brightness  value  arrays  that  were 
centered  on  stations  with  large  temporal  variations  in  snowfall 
intensity  was  formed.  Four  snowstorm  cases  were  selected  from  a  t~tal 
of  fifteen  cases  examined  during  the  winters  of  1979-80  and  1980-81. 
Selection  was  based  upon  the  criteria  of  variability  of  snowfall 
intensity  and  cloud  structure  during  the  storms.  A  total  of  210  arrays 
were  used,  each  having  a  required  standard  deviation  in  brightness 
values  of  at  least  2.0, 

Five  numerical  parameters  were  then  estimated  from  each  array  of 
digital  IR  brightness  values.  They  included  deviation  of  station 
brightness  (DM)  from  the  average  array  brightness,  two  five-point 
station  centered  Laplacian  of  the  brightness  field  (one  with  a  16  nm 
horizontal  scale  (ILAP1)  and  the  other  with  a  32  nm  horizontal  scale 
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(0LAP2 ) )  and  two  nine-point  station-centered  Laplacians  of  the 
brightness  values,  with  horizontal  scales  the  same  as  the  five-point 
structure. 

Mean  values  and  interval  frequency  histograms  were  constructed  for 
each  of  the  five  parameters  according  to  observed  snowfall  intensity  at 
the  array  centered  station.  Results  based  upon  the  distributions  of 
the  means  and  frequencies  showed  that  for  the  most  skillful  indicator 
of  intensity  (0LAP1),  positive  determination  of  precipitation  intensity 
at  a  locale  was  possible  only  approximately  10  per  cent  of  the  time. 

The  hourly  changes  in  the  five  parameters  was  examined  in  a  similar 
manner  for  cases  of  intensity  (visibility)  change  from  hour  to  hour. 

All  parameter  changes  displayed  distributions  expected  from  an 
Idealized  relationship  of  cloud  structure  to  intensity.  The  most 
skillful  indicator  of  intensity  changes  (aOLAP2)  was  able  to  correctly 
infer  the  change  nearly  30  per  cent  of  the  time. 

Based  upon  these  results,  inference  of  intensity  for  cases  of 
continuous  precipitation  in  middle  latitudes  does  not  seem  feasible  if 
only  relative  variations  in  the  IR  brightness  values  distributed  over  a 
given  location  are  considered.  Intensity  change,  however,  apparently 
can  be  inferred  a  significant  percentage  of  the  time,  and  study  aimed 
at  further  developing  these  techniques  seems  appropriate. 

Several  problems  with  the  methods  used  were  realized  during  the 
course  of  the  study.  The  primary  problem  seemed  to  be  horizontal 
displacement  of  snow  during  its  fall  from  its  level  of  origin  to  its 
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point  of  deposition;  a  problem  that  could  not  be  compensated  for  in 
this  research.  The  observed  Intensity  at  the  surface  and  the  location 
of  the  clouds  producing  the  heavier  precipitation  in  these  cases  of 
strong  wind  shear  were  not  colocated.  In  a  preliminary  comparison, 
using  Case  1  data,  maxima  fields  of  the  radar  return  and  satellite  IR 
brightness  were  co-located  10  to  15  miles  northeast  of  the  surface 
points  of  maximum  snowfall.  If  the  model  had  included  provisions  for 
the  drift  of  falling  snow  due  to  wind  sho.r,  the  evidence  suggests  that 
the  discrimination  would  have  been  increased  substantially.  In  future 
studies  that  address  the  nowcasting  problem  of  mesoscale  precipitation 
in  the  form  of  snow,  a  degree  of  freedom  should  be  included  to  account 
for  the  displacement  of  falling  snow  of  convective  origin.  This 
problem  would  undoubtedly  be  much  less  prominent  for  rainfall  cases. 

A  technique  which  centers  the  arrays  on  easily  tracked  brightness 
maxima  and  minima  is  suggested.  The  arrays  would  be  moved  to  coincide 
with  positions  of  the  maxima  or  minima  and  the  discrimination  ability 
of  parameter  values  and  parameter  changes  would  be  tested  for  stations 
located  near  and  downwind  from  the  array  centers. 

Another  problem  .vas  the  inability  of  the  conventional  five-point 
l.aplacian  structure  to  adequately  describe  variations  in  the  IR 
brightness  field  for  banded  convection.  This  problem  was  reduced  by 
using  a  nine-point  Laplacian  structure  (eight  perimeter  points  and  one 
central  point). 

The  same  methods  to  estimate  snowfall  intensity  change  were 
applied  to  a  special  array  subset  of  fifty  occurrences  of  increased 
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variability  of  the  brightness  values  within  the  arrays.  The 
reliability  of  estimation  of  snowfall  Intensity  changes  from  the  five 
parameter  was  not  increased  In  this  subset. 

In  retrospect,  the  use  of  only  snowstorm  cases  in  this  research 
was  likely  a  limitation  In  the  inference  of  precipitation  intensity 
from  satellite  information.  Variations  in  observed  intensity  of 
snowfall  are  possibly  not  as  well  related  to  parameters  used  in  this 
study  as  the  relation  might  be  with  rainfall.  Rainshowers  in  warmer 
seasons  are  associated  with  a  warmer  and  somewhat  less  stable 
troposphere.  Thus,  depth  and  organization  of  the  convective  systems  of 
rainshowers  is  likely  more  pronounced.  Consequently,  parameters 
utilizing  cloud  top  temperatures  to  estimate  showery  precipitation 
should  discriminate  better  in  conditions  of  rainshowers  than 
snowshowers.  In  addition,  the  problem  of  displacement  would  be  reduced 
due  to  the  increased  fall  velocity  of  raindrops  over  snowflakes. 
However,  because  of  these  different  fall  rates,  the  analysis  to  infer 
precipitation  intensity  and/or  intensity  changes  should  be  completed 
separately  for  snow  and  rain,  due  to  the  ambiguity  of  larger 
displacement  for  snow  in  a  wind-sheared  atmosphere. 
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gure  1:  Typical  distributions  of  mesoscale  convective  bands 
associated  with  (a)  an  occluded  cyclone  and  (b)  a 
developing  cyclone.  Arrows  indicate  movement  of 
convective  activity. 


Figure  2:  NMC  surface  analyses  for 
(a)  1200  GMT  14  April  1980,  (b) 

1800  GMT  14  April  1980  and 

(c)  0000  GMT  15  April  i960.  Stars 
in  (a)  locate  stations  on  which 
brightness  arrays  were  centered, 
circles  indicate  sounding  stations 
used. 
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Figure  3:  NMC  500  mb  geopotential  analyses  (heights  in  decameters)  for 
(a)  1200  GMT  14  April  1980  and  (b)  0000  GMT  15  April  1980. 
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Figure  4:  Vertical  distributions  of  temperature  (solid  lines)  and 
dewpoint  (dashed  lines)  for  (a)  Peoria,  Illinois  and  (b) 
Green  Bay,  Wisconsin,  1200  GMT  14  April  1980. 
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8:  Same  as  Figure  3  except  for  1200  GMT  6  February 
1980. 


9:  Same  as  Figure  4  except  for  (a)  Dulles  International 
Airport,  Maryland  and  (b)  Greensboro,  North  Carolina 
1200  GMT  6  February  1980. 
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Figure  13:  Same  as  Figure  3  except  for  (a)  1200  GMT  2  March  1980 

and  (b)  0000  GMT  3  March  1980. 


Figure  14:  Same  as  Figure  4  except  for  (a)  Greensboro,  North 

Carolina  and  (b)  Cape  Hatteras,  North  Carolina,  1200 
GMT  2  March  1980. 


Figure  18:  Same  as  Figure  3  except  for  (a)  1200  GMT  31  January  1981 
and  (b)  0000  GMT  1  February  1981. 
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Figure  19:  Same  as  Figure  4  except  for  Omaha,  Nebraska,  1200 
GMT  31  January  1981. 
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Figure  20: 


Example  McIDAS  listings  of  weather  observations,  6 
February  1980  for  Raleigh-Durham,  North  Carolina. 
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Fiqure  21:  (a)  Sample  digital  IR  brightness  array  (centered  station  is 
at  x)  with  (b)  mean,  median,  mode  and  standard  deviation 

of  values. 
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Figure  22:  Example  of  point  brightness  calculation  (for  station 
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Figure  23:  Example  showing  magnitude  dependency  of  Laplacian  values 
upon  the  "circularity"  of  relative  maxima  or  minima. 
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figure  24:  Example  showing  dependency  of  Laplacian  values  upon 
the  orientation  of  banded  maxima  or  minima. 
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LAP  2.  —  B^ + B2+  B3+  B4+B5 
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Figure  25:  Schematic  of  five-point  laplacians  (ILAP1  and  0LAP1) 
and  nine-point  Laplacians  (ILAP2  and  0LAP2) 
including  equations  used  for  calculations  (x  and  y  are 
16  nm  or  32  nm). 


Figure  26:  Tabulation  format  for  hourly  values  of  visibility, 
weather  and  all  five  parameters  (sample  data  for 
Madison,  Wisconsin  from  14  April  1980  case). 
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Figure  27: 


Tabulation  format  for  hourly  changes  in  visibility 
and  all  five  parameters  (example  for  intensity 
increases). 
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Figure  33:  Same  as  Figure  32  except  for  (a,  b)  AllAPl  and  (c,  d)  aOLAPI 
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Figure  43:  Same  as  Figure  42  except  for  (a)  A11AP1,  (b)  aOLAPI,  (c) 
^IL4P2  and  ( d )  A01AP2. 


